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A general graded reflection equation algebra is proposed and the corresponding boundary quan- 
tum inverse scattering method is formulated. The formalism is applicable to all boundary lattice 
systems where an invertible R-matrix exists. As an application, the integrable open-boundary con- 
ditions for the g-deformed supersymmetric U model of strongly correlated electrons are investigated. 
The diagonal boundary K-matrices are found and a class of integrable boundary terms are deter- 
mined. The boundary system is solved by means of the coordinate space Bethe ansatz technique 
and the Bethe ansatz equations are derived. As a sideline, it is shown that all R-matrices associated 
with a quantum afflne superalgebra enjoy the crossing-unitarity property. 
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In the last decade, much attention has been paid to the study of strongly correlated electron systems. In particular, 
the discovery of high-T c superconductivity has greatly stimulated investigations on various electron lattice models in 
one dimension (ID), which are exactly solvable by means of the coordinate Bethe ansatz method or the quantum 
inverse scattering method (QISM). The most known and studied integrable electron systems are perhaps the Hubbard 
£^ . and the supersymmetric t-J models. Other integrable correlated electron theories of interest include, for instance, the 
so-called extended Hubbard model (0 , the Hubbard- like model M , the supersymmetric U model || and its eight-state 
version and the g-deformed supersymmetric U model H§|. These models have been extensively investigated in 
the literature (c.f.: 

fN| , On the other hand, one of the recent developments in the theory of completely integrable lattice models is Sklyanin's 
work [[[o| on the boundary QISM, which may be used to describe integrable systems on a finite interval with inde- 
pendent boundary conditions on each end. The important ingredient in this boundary QISM is the new algebraic 
structure, the reflection equation (RE) algebra. The solutions to the REs are called boundary K-matrices which in 
turn give rise to boundary conditions compatible with the bulk model integrability. 

In Sklyanin's formulation on the boundary QISM, the R-matrices are assumed to enjoy P-, T- and crossing- 
£ — symmetry. Sklyanin's results are generalized by Mezincescu et al pr| to treat cases where, instead of having the 
OS . separate P- and T-symmetry, the R-matrices satisfy, in addition to the crossing-symmetry, the less restrictive condition 
of the combined PT invariance. It is noted by de Vega et al jl^] that the formalism of Mezincescu et al is actually also 
applicable even if the R-matrices only have the so-called crossing-unitarity, a weaker version of the crossing-symmetry. 

Generalizations to the following two cases seem very interesting: (i) the graded or supersymmetric case, and (ii) 
the case where the R-matrices do not obey any constraint conditions (except the unitarity condition). Many attempts 
have been made in the literature concerning the extension (i). However, very few authors treat the grading properly 
from the beginning to the end. Therefore, in the eyes of the present authors most known results in the literature are 
^ not very satisfactory and a fully graded or supersymmetric formalism is desirable. Concerning (ii), two important 
examples are the R-matrices [pi corresponding to the models of two coupled and three coupled XY spin chains 
. £h ! introduced in , to which the formalism developed in does not apply. 

In section II of this paper, we shall fulfill the two aims of extensions. More specifically, we shall formulate a fully 
supersymmetric boundary QISM which is applicable to any cases where an invertible R-matrix exists. We introduce 
a very general graded RE algebra and show that this algebra indeed leads to a commuting family of the boundary 
transfer matrices. Throughout the procedure, no spectral parameter multiplicativity/additivity of the graded quantum 
Yang-Baxter equation (QYBE) has been assumed and no constraint conditions on the R-matrices (except the always 
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satisfied unitarity property) been imposed. Our formalism is a supersymmetric generalization of that developed in 
|l3].(i4| for the bosonic (or non-supersymmetric) case, where open-boundary integrable conditions for the models of 
two coupled and three coupled one-dimensional XY spin chains have been constructed, respectively. We then use, 
in section III, our results to study integrable open-boundary conditions for the g-deformed supersymmetric U model 
|!|,[| of strongly correlated electrons. The quantum integrability of the model in the bulk has been established in || 
by embedding the bulk model Hamiltonian into a one-parameter family of commuting transfer matices formed from 
a E/q[<?Z(2|l)] invariant R-matrix. We solve the graded REs for the diagonal boundary K-matrices and determine the 
open-boundary integrable model. The boundary model Hamiltonian is shown to be related to the second derivative 
of the boundary transfer matrix. In section IV, we solve this boundary model by the coordinate space Bethe ansatz 
method and derive the Bethe ansatz equations. In the Appendix, we show that all R-matrices associated with a 
quantum affinc superalgebra possess the crossing-unitarity. 



II. GRADED REFLECTION EQUATIONS AND TRANSFER MATRIX: GENERAL FORMULATION 

In this section, we establish a very general RE algebra. We shall not assume spectral-parameter multiplicativ- 
ity/additivity of the graded QYBE, although for the special example we shall discuss in section III the spectral 
parameters are multiplicative/additive. Throughtout the procedure, no constraint conditions have been assumed on 
the R-matrices. 

To begin with, let V be a finite-dimensional Z2 graded linear superspace. Let R G End(V ® V) be a solution to 
the graded QYBE 

Rl2(ui,U2)Rl3(u 1 ,U 3 )R 2 3(u2,U3) = i?23(«2, U 3 )R 13 (u 1 ,U 3 )Ri 2 {u 1 ,U 2 ). (II. 1) 

Here Rjk{u) denotes the matrix on V (g) V <g> V acting on the j-th and fc-th superspaces and as an identity on the 
remaining superspace. The variables u\, u 2 and U3 are spectral parameters. The tensor product should be understood 
in the graded sense, that is the multiplication rule for any homogeneous elements x, y, x' ', y' G V is given by 

(x ® y)(x' <8> y') = (-l) lslM (xx 1 ® yy') {11.2) 

where [x] stands for the grading of element x: i.e. [x] = if x is even (or bosonic) and [x] = 1 if a: is odd (or fermionic). 
Let P be the graded permutation operator in V ® V. Then P(x ® y) = (—l)™ y ^y <8> x for all homogeneous x, y G V 
and i?2i(u) = P\2Ri2{u)P\2- We form the monodromy matrix T(u) for an L-site lattice chain 

T(u)=L QL (u)---L Q1 (u), (II.3) 

where Loj(u) = Roj(u,0), where subscript labels the auxiliary superspace V. Indeed, one may show that T(u) 
generates a representation of the graded quantum Yang-Baxter algebra, 

Ri2{ui,u 2 ) T (111) T (u 2 ) =T (u 2 ) T (u 1 )R 12 {u 1 ,u 2 ), (II.4) 

1 2 
where X= X <8> 1 and X= 1 <8> X, for any matrix X G EndiV). 

In order to construct integrable electronic models with open boundary conditions, we need to introduce an appropi- 

ate graded RE algebra. We introduce two associative superalgebras 71 and T + defined by the R-matrices and the 

relations 

1 22 1 

Ru(ui,U2)T- (ui)R2i(u2, -ui) T- {u 2 ) = T- (1*2)^12(^1, -«2j T- {u 1 )R 2 i(-U2,-u 1 ) 7 

R$ ^{u2i Ul ) Tf 1 (u 1 )R 12 (-u 1 ,u 2 ) Tl st - (ua) = (ua^if-ua.ui) Tf 1 {m)Rf 2 l lst2 (- Ul ,-u 2 ), (11.5) 
where we have defined new objects R and R through the relations 

Rf 2 2 {-Ul,U2)R'2i {u 2 , -Ml) = 1, 
■zisti . 

R 2X (-U2,«i)i^ a (ui,-u 2 ) = l, (II.6) 

and, as defined in the Appendix, sti stands for the supertransposition taken in the i-th space, whereas isU is the 
operation inverse to sti. As will become clear below, one of the important steps towards formulating a correct 
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formalism for the graded or supersymmetric case is to introduce, in the second RE in (II. 5) below, the inverse 
operation of the supertransposition. The introduction of this inverse operation is essential because, applying the same 
supertransposition operation twice does not in general give an identity operation. For instance, applying si,, i = 1, 2 
, twice to the R-matrix R(u) does not yield R(u): 

{R{u) st >} st > + R(u), i = 1,2. (11,7) 

In all cases, the quantum R-matrices possess the unitarity [], 

i?i 2 (ui,u 2 )i?2i(u2,ui) = 1- (H.8) 



We now show that the second RE in (II. 5) is indeed the correct "conjugation" to the first one, so that the boundary 
transfer matrices defined as usual constitute a commuting familty. Following Sklyanin's arguments fLOfl , one defines 
the boundary transfer matrix t(u) as 

t(u) = str(T+(u)T-(u)), {11.9) 

where str denotes the supertrace taken over the auxiliary superspace V. Then it can be shown that 

[t(«i),t(u 3 )]=0. (11.10) 

The proof is elementary. We nevertheless present the details. By means of the commutativity of operators in 7+ and 
71, 

11 2 2 

ii(ui)i 2 (u 2 ) = stri{T + (tti) 71 (u 1 )}str 2 {T + (u 2 ) T_ (u 2 )} 

- st n {T^ ( Ul ) T_ stl (u 1 )}str 2 {T + (u a ) T_ (u 2 )} 
= str 12 {Tf^ ( Ul ) T+ (7i 2 ) Ti h (ui) T_ (u 2 )} = ■■ ■ . 



Inserting a variant of the first expression in ([1.6) for "1" into the supertrace 



l 2 sf i 2 

• • • = strx 2 {T+ 1 (ui) T + (u2)i?i2 2 (-ui,u2)^2i 1 ( u 2, -ui) T"- 1 (ui) T_ (u 2 )} = • • • , 

then applying the supertransposition, 

••• =str 12 {Tf l {u x )R l2 {-ui,u 2 )Ti st * M} st2 {T- [u x )R 21 {u 2 ,-u x )T^ (u 2 )} stl 
= str 12 {T* h {u x )R l2 {- Ul ,u 2 ) Tl st * (u 2 )} istl S ' 2 {T- (u x )ii2i(« 2 ,-«i) T_ (u 2 )} 

" * * 5 

and, finally, inserting the unitarity property into the supertrace one obtains 

•■• =str l2 [{T^ ( Ul )R 12 (-ux,u 2 )Tl st * (u 2 )} 1 ^ st > R 21 (u 2lUl ) 

1 2 

x fli 2 (ui,« 2 ){71 (ui)-R 2 i(u 2 ,-Ui) T_ (u 2 )}] 
= str 12 [{R% ist2 (u 2 , Ul ) Tf 1 (uiJ.RiaC-ui.Ua) T^ 2 M}"' 1 s * 2 

1 2 

x J?i 2 (mi,« 2 ) T_ R 2 i(u 2 , -ui)(ui) T_ (u 2 )] = ••• . 



Applying the RE algebra ( |tl.5|) , 



1 One can always normalize the R-matrices so that the right hand side of (H.8) equals to 1 
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1 

-sty 



str 12 [{Tl st > (u 2 )R 21 (-u 2 , Ul ) T + Stl ( Ul )R s £ ist2 (- Ul ,-u 2 )} ist i st * 

2 1 

x T_ (u 2 )Ri2(ux,-u 2 ) 71 {ux)R 21 {-u 2 , -m)] 
str 12 [Tl st2 {u 2 )R 2l {-u 2 ,u x )T^ ( Ul )R s £ lst H~ Ul , -u 2 ) 
x {71 (u 2 )iii2(ui,-U2)r_ (uiJiZaif-ua.-ui)}** 1 ls * 2 ]. 



Repeating the whole chain of transformation in the reverse order and keeping in mind the second expression in (EI. 6) 
for "1", one ends up with t(u 2 )t(ui), as required. 

One can obtain a class of realizations of the superalgebras 7+ and 71 by choosing T± (u) to be the form 

T-{u) = T_(u).K'_(u)T:: 1 (-u), T+{u) = K+{u) 



(11.11) 



with 



T_(u) = R 0N (u) ■ ■ ■ Roi(u), 



(11.12) 

where iV is any node between 1 and L (including 1 and L), and K±(u), called boundary K-matrices, are representations 
of T± in a Grassmann algebra. In this realization, the elements of the matrices K±(u) are all Grassmann numbers. 
Note that the K-matrices K±(u) sati sfy the same relations as T±(u), respectively. That is, the K-matrices fulfill the 
following graded REs deduced from flll.5| ) 

1 2 2 1 

Ri 2 {ui,u 2 ) K- (ui)R 2 i(u 2 , -mi) K- (u 2 ) = K- {u 2 )R\ 2 {u\, -u 2 ) K- (u 1 )R 2 i(-u 2 , 



R% ™ t2 (u 2 , Ul ) K s ^ ( Ul )R 12 (- Ul ,u 2 ) K'f 2 («a) = K'f 2 (u 2 )R 21 (-u 2 , Ul ) K s f ( Ul )R% lst2 -u 2 ). (11.13) 



The REs ( P-5|) and their realizations ( |II.13|) in the Grassmann algebra are generalizations of those introduced in 
|^3|,^6[ . Note that no constraint conditions have been imposed on the R-matrices. Therefore our graded REs apply 
to any case where an invertible R-matrix exists. 

Let us now consider the important "special" case in which the R-matrices are related to finite dimensional represen- 
tations 7ry of the quantum affine superalgebra U q [Q^ (k = 1, 2) for generic q, where Q is any simple Lie superalgebra. 
In the Appendix, we show that all such R-matrices enjoy the crossing-unitarity properties [see (All) or (A.14)]. One 
can show that there exists a more general realization of the superalgebra T + 



T + (u) = R 0L (u) ■ ■ ■ R ,n+i (u) ■ 



(11.14) 



In the sequel, without loss of generality, we shall choose N = L so that T + (u) = K + (u). As the spectral parameters 
in the QYBE are now multiplicative/additivity, the K-matrices fulfill the REs of the following form, 

1 2 2 1 

R\ 2 {ui - u 2 ) K- (ui)i?2i(«i + u 2 ) K- (u 2 ) = K- (u 2 )Ri 2 (ui + u 2 ) K- {ui)R 2 i{u 1 - u 2 ), 
R 12 {- Ul +u 2 ) K+ ( Ul )R 2 T St2 (- Ul -u 2 ) K+(u 2 ) = K+ {u 2 )R^ St2 {- Ul - u 2 ) K+ {u x )R 21 {-u x + u 2 ), (11.15) 
where in deriving the second relation we have applied the operation ist\ st 2 on both side of the second equation of 

dp . 



By (II.6) 



R 2 T S * 2 (-^i - u 2 ) = (((R 21 (- Ul - u^ 1 )^)- 1 )^- 
(((iZuC-ui-t^)- 1 )** 1 )- 1 )"* 1 



fiisti St2 

/l 12 



(-ttl 



U 2 



(11.16) 



With the help of (A. 11), (A. 9) and (A. 12) (one should identify z = q u t o convert the multiplicative spectral parameter 
z into the additive one u), one can show that the second RE in ( 1 1 . 1 5| ) becomes 



l 1 2i 2 

iM-«i + "2) K+ ( Ul ) A!' 1 R 2X {-ui -u 2 + -g) MK+ (u 2 ) 



1 2 

--MK+ (u 2 )R 12 ( 



2 i l 
-U\ — u 2 + —g) M K + (u 1 )R 21 (—u 1 



u 2 ), 



(11.17) 



where M = irv{q~ 2p ) is the so-called crossing matrix and g is defined as in the Appendix. The RE (11.17) coincides, 
in the case of k = 1, with the one used in and is applicable to all cases whose R-matrices are related to the 
quantum affine superalgebra U q [Q^ k '] for generic q. 
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III. INTEGRABLE BOUNDARY K-MATRICES AND Q-DEFORMED SUPERSYMMETRIC U MODEL 

WITH BOUNDARY TERMS 



Let Cj, a and c t a denote fermionic creation and annihilation operators with spin a at site j, which satisfy the 

anti-commutation relations given by {c\ a ,Cj^ T } — <5y<5 CTT , where i,j — 1,2, ■■■,L and er, r =f, |. We consider the 
g-deformed supersymmctric U model with boundary terms described by the following Hamiltonian: 

L-l 

Bl + B R , 



H = J2 H 



3,3+1 



(III.1) 



3=1 



where Hjj + i is the bulk Hamiltonian density of the g-deformed supersymmetric U model 



3+1, 



+ -^\ n 3\ n 3l+ n 3 + l\ n 3 + l\) 

+t p( c n c n c 3+n c 3+n + h - c -) + f 1 n j + M -1 n j+1 , 

and Bl and Br are left and right boundary terms, respectively, given by 

/i — /i -1 / sinh7 



(III.2) 



Br 



Br 



2 sinh 



7(2~g-) 



sinh 



-ni T n u 



2 

sinh 7 



2 sinh 



7(2~g+) 
2 



sinh 



7£+ 



where n J(J is the density operator n, CT = c 



ja 3 a , 



% 3l 



n L i i n Ll - e 



riji and 



-7(1- 



7 (1- 



'n L 



(III.3) 



t p = — = e[2e ^ (cosh 77 — cosh 7)] 2 , e 



sinh^ — 7)/2 



sinh(r? + 7)/2' 



(III.4) 



£± are some parameters describing the boundary effects. It is interesting to note that the Hamiltonian (III. 1 ) becomes 
[/ 9 [<?^(2|l)]-invariant in the limits — ► —00, — ► 00. 

We shall establish quantum integrability for the system defined by the Hamiltonian (111.1), by applying the general 
formalism developed in the previous section. To this end, let us first of all recall some basic results of the g-deformed 
s uper symmetric U model with the periodic boundary conditions. In it was shown that the bulk Hamiltonian 
(III. 2) of the model commutes with the bulk transfer matrix r(u), which is the supertrace of the monodromy matrix 
T(u) with the local monodromy matrix L$j(u) — Roj(u). That is, 



T (u) = str(T(u)), T(u) = R QL {u) ■ ■■R 01 (u), 
where the quantum R-matrix R(u) = PR(u) 1 with 



R{u) 



„u _ 2a 1 _ „u+2a+2 

' ' -P1+P2+ - — P 3 , 



1 — q u + 2a 



qu _ g2a+2 



(111.5) 



(III.6) 



where Pi, i = 1, 2,3 are the projection operators whose explicit formulae may be found in 0], 

In order to describe integrable systems with the boundary conditions different from the periodic ones, we first solve 
the REs for the two boundary K-matrices K±(u). For our purpose, we only look for solutions where K±(u) are 
diagonal. After complicated algebraic manipulations, we find 



K-{u) 



sinh sinh 7 ^ 2 2 - ) 



/A_(u) 

B_(it) 

B_{u) 

\ C_(u) 



(HI.7) 
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where 



A-{u) 



' sinh 



^- + ^sinh^- 2 + ^. 



w 2 2 

C_ («) = e-T« sinh ~ U) sinh 7(-"~2 + g-) 



(III.S 



As is shown in the Appendix, the R-matrix (III. 6) satisfies the crossing-uniatriy condition. This implies that there is 
an isomorphism between the graded REs for K + and K- : 



K+(u) = MK-(-u + l), 
where M is given by (up to an overall factor) 



M = 



10 

10 

e 2 ^ 

e 2 ^ 



(III.9) 



(111.10) 



Therefore we may choose the boundary K- matrix K + (it) as 

K+(u) = 



A + {u) 

B+(u) 

C+(u) 






(111.11) 



with 



e-rC-D sinh 7(2 - 2a - £ + - u) 7 (2a + £ + + „) 
2 2 

-27 • y, 7(-2«-C++^) . , 7(2« + g+ + u) 
e ' sinn Sinn , 



A + (u) 
B+(u) 

C + (u) = sinh 7(-2q-g + +«) Binh 7(2a + g + A) 



U + («) = sinh 7(~2a sinh 



7(2 + 2a + £+ - u) 



(III. 12) 



Indeed,we have checked that this K + matrix constitutes a solution to the graded REs JII.17 ). 

To show that the Hamiltonian ( III. 1 ) can be embedded into the boundary transfer matrix t(u) constructed in section 
II is an involved algebraic manipulation. This is because the supertrace of K + (0) is equal to zero. So at best we 
can only expect that the Hamiltonian (III.l) appears as the second derivative of the boundary transfer matrix with 
respect to the spectral parameter u, at u = 0. 

Let us expand the local monodromy matrix LqAu) up to the second order in the spectral parameter u, 



Lq 3 {u) = (1 + H j0 u + yBjou 2 + + ■■ -)L 0j (0). 



(III. 13) 



Substituting this expression into the boundary transfer matrix t(u), and after a lengthy but straightforward algebraic 
calculation, one finds 



t{u) = C x u + C 2 (H + const.)u 2 



(III. 14) 



where Cj(i = 1, 2, • • •) are some scalar functions of the boundary constant £+. Then it can be shown that up to some 
additive constants the Hamiltonian (tll.l) is related to the second derivative of the boundary transfer matrix, 



H 



q a+1 - i~ a '\ H R 



In q 



t"(0) 



4(V + 2W) 



2(V + 2W) 



str K+ (O)Glo 



+2 str Q K'+ (0)fff o + str K+ (0) (H* 0) 



(111.15) 



G 



where 



V = str K' + (0), W = str (#+ (O)if£ 



(111.16) 



if we make the following identifications: 



7 Q+1 



<r 



(111.17) 



Thus, we have shown that the Hamiltonian (1II.1) of the g-deformed supersymmetric U model with the boundary 
terms Bl and Br is related to a class of commuting transfer matrices. As a result, the system has an infinite number 
of higher conserved currents which are involutive with each other, and therefore the system under study is completely 
integrable. 



IV. THE BETHE ANSATZ EQUATIONS 

Having established the quantum integrability of the model, let us in this section solve it by using the coordinate 
space Bethe ansatz method. Following |17| , ^6| , pt , we assume that the eigenfunction of Hamiltonian (III.l) takes the 
form 

l*}= ^^,-,a N (xi,---,X N )cl ia . l ---cl NaN \0), 

N 

*ffi,-,<7jv (xx, ■ ■ ■ , x N ) = Y e P A (T Q i,-,a QN {kpQi,- ■ ■ i k PQN) exp(i 2J kpjXj), (IV.l) 

P 3 = 1 

where the summation is taken over all permutations and negations of kx, ■ ■ • , A; at, and Q is the permutation of the N 
particles such that 1 < xq% < • • • < xqn < L. The symbol ep is a sign factor ±1 and changes its sign under each 
'mutation'. Substituting the wavefunction into the eigenvalue equation H\^>) = E\^f), one gets 

A... f aj ,&%,••• (' ' ' 5 kj , • ■ ■) — Sij {ki , kj ) A. . . ai j(7 . ; . . . (• • • , k i , kj , • • •) , 

kj , • • •) = S L (kj\ pi ai )A l j i ^...(kj, ■ ■ •), 

A.., CT< (- • • , -kj) = s R (k f ,p L(Ti )A... tai (- ■ • , kj), (IV.2) 
where Sij(ki,kj) are the two-particle scattering matrices, 

Sij(ki,kj)n = Sij(ki, kj)<22 = 1, 

1 (k k ) 12 - <?•(*•• k) 21 - sm ( A » ~ Aj) 

- fcj ai - gin(Ai _ A i<y j. 



Sij[ki, kj) 21 



j 

12 p _ i(Ai _ A .) S1TU7 

sin(Ai — Aj — «7) ' 

with Aj being suitable particle rapidities related to the quasi-momenta kj of the electrons by |5| 

0(A,a) = 2arctan(cothatanA), 

1 f r sinhi(77 + 7) n 1 
° = 7 ln[ i ] - 7 , (IV.4) 



sinh ^ (?7 — 7) 

and s L (kj-,pi ai ), s R (kj;pL ai ) are the boundary scatering matrices, 
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s L (kj]p lai ) 



s R {kj-,p Lo 



1 -y\ ai e iki 
1 -pi ai e~ lkj 

1 -Vhofi %k i 



(IV.5) 



with 



Via = Pi 



PLa = PL 



-7(1"^) 



2 sinh 



7(g-|-) 



2 sinh 



t(2-|+) 



=7(1-^) 



(IV.6) 



Then, the diagonalization of Hamiltonian (tll.l) reduces to solving the following matrix eigenvalue equation 

Tjt = t, j = l,---,N, 
where t denotes an eigenvector on the space of the spin variables and Tj takes the form 

Tj = Sj (k^i-k, ; Plaj )RJ (k J )R+ (kj)s R (kj ; p LrTj )S+ (k, ) 

with 

= Sj.x{kj, k N ) ■ ■ ■ Sjj + i(kj, fcj+i), 
Sj (%) = Sjj-i(kj, kj-i) ■ ■ ■ Sj ! i(kj, fci), 
Rj(kj) = Si,j(ki,-kj) ■ ■ ■ Sj-ij(kj-i,-kj), 
Rffc) = S j+ i,j{k j+ u ~kj) ■ ■ ■ S Ntj {k N , -kj). 

This problem can be solved using the algebraic Bethe ansatz method. The Bethe ansatz equations are 

^ 2(i+1) C(^;Pi)C(^L) = ~ {Xj A ° + f } Sin(Aj + Aq + f } 



^ sin(Aj - A Q - §) sin(Aj + A Q - f ) 

^ sin(A Q - Xj + ^) sin(A Q + Xj + ^-) _ yr sin(A Q - Ap + iy) sin(A Q + A^ + by) 
4i sin(A Q - A~ - f ) sin(A Q + A~ - f ) ~ ji, sin(A Q - A/3 - i 7 ) sin(A Q + A - i 7 ) ' 



(3 = 1 



(IV.7) 



(IV.8) 



(IV.9) 



(IV.10) 



where ((k;p) = (1 — pe lk )/(l — pe lk ). The energy eigenvalue E 1 of the model is given by E = — 2^. =1 cosfcj (up to 
an unimportant additive constant, which we have dropped). 



V. CONCLUSION 



We have proposed a very general graded RE algebra and developed the corresponding boundary QISM. In formu- 
lating our general formalism in section II we do not impose any constraint conditions on the quantum R-matrices, and 
therefore our formalism applies to all lattice boundary systems where an invertible R-matrix exists. Two nontrivial 
examples are the supersymmetric versions of the models of two coupled and three coupled one-dimensional XY spin 
chains Ml, where the fermionic R-matrices do not possess the so-called crossing- unitarity and so no isomorphism 
between the boundary K-matrices K- and K + exists. We have also considered the important "special" case in which 
the R-matrices are associated with finite-dimensional representations of the quantum affine superalgebra U q [G^] for 
generic q. For such a case the two REs are isomorphic to each other since, as is shown in the Appendix, all such 
R-matrices enjoy the crossing-unitarity property. 

We have applied our general formalism to study integrable open-boundary conditions for the g-deformed supersym- 
metric U model of strongly correlated electrons. The quantum integrability of the boundary system is established by 
the fact that the corresponding Hamiltonian may be embedded into a one-parameter family of commuting transfer 
matrices. Moreover, the Bethe ansatz equations are derived by means of the coordinate space Bethe ansatz approach. 
This provides us with a basis for computing the finite-size corrections to the low-lying energies in the system, which 
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in turn allow us to use the boundary conformal field theory technique to study the critical properties of the boundary 

model. 

As mentioned in section III, the Hamiltonian (HI.l) becomes {7 g [gZ(2|l)]-invariant in the limits £_ — > — oo and 
£ + — > oo. It is very interesting to see whether or not the Bethe states constructed above constitute the highest weight 
states for t/ g [gZ(2|l)]. This property is crucial in understanding the completeness of the Bethe states. A similar 
problem has been studied |L8| for the J7 g [<7Z(2|l)]-invariant open t — J chain. It seems interesting to note that although 
the completeness problem for the Bethe states are well studied in the periodic case for the Hubbard model and the 
supersymmetric t — J model fl9[| , it remains largely unexplored in the non-periodic (or open) boundary case. Another 
interesting question is to derive the Bethe ansatz equations using the algebraic Bethe ansatz approach. 
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Grant and External Support Enabling Grant. H.-Q.Z would like to thank the Department of Mathematics, University 
of Queensland, for kind hospitality. He is also partially supported by the National Natural Science Foundation of 
China and Sichuan Young Investigators Science and Technology Fund. 

APPENDIX A: ON THE CROSSING-UNITARITY 

So far in literature, the discussions on the crossing-unitarity of R-matrices have been on a case by case basis and 
the crossing-unitarity of a given R-matrix has been checked by brute force. 

In this appendix, we show that all R-matrices associated with finite-dimensional representations of the quantum 
affine superalgebra U q [G^] [k — 1, 2) for generic q, where Q is any simple Lie supcralgcbra, enjoy the crossing unitarity 
property. 

Let us first of all recall some facts about the affine superalgebra Q^ k \ Let Go be the fixed point subalgebra under 
the diagram automorphism f of Q of order k. In the case of k = 1, we have Go = Q. For k = 2 we may decompose G 
as Go plus a ^-representation Gi of G- Let 

f highest root of Go = G, for k = 1, 

1 highest weight of the Go ~ representation Gi, for k = 2. \ ■ ) 

Following the usual convention, we denote the weig ht of G [k) by A = (A, k,t), where A is a weight of Go ■ With this 
convention the nondegenerate form ( , ) induced on the weights can be expressed as 

(A, A') = (A, A') + kt + k't. (A.2) 

Let p denote the distinguished dominant weight of G^ satisfying (p, ai) = \{oti, cxi), where a, (0 < i < r) are simple 
roots of G^ ■ Then p is given by 

P = P + 91, 7 =(0,1,0), (A.3) 

where g = f ("0, ip + 2p) and p is the graded half sum of positive roots of Go- 

We shall not give the defining relations for U q [G^], but mention that the actions of coproduct and antipode on its 
generators {hi, e^, /j, d, < i < r} are given by 

A(hi) = hi<8>l + l®hi, A(d) = + 

A(ef) = e t O +q^ ® e h A(/ 4 ) = / t <g> q^ + q^ ® fa 

S(a) = -q-Pa <f , Va = d, h, e h fa (A.4) 
Define an automorphism D z of U q [G^] by 

D z {ei) =z 5 *°e u D z {fi)=z- Si0 fi, D z {K) = h h D z {d) = d. (A.5) 

Then 

S 2 (a)=q- 2 "D q _ ig {a) q 2 ", ^ 2 (a) = q 2 " D qh (a) q~ 2 ", Va e U q {G { % (A.6) 

which can be checked on the generators [remembering that the simple roots associated with eo, fo are ao — ±(r5— if)), 
respectively, where S = (0, 0, 1)]. We define the right dual module V* and left dual module *V of V by 
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7r y . (a) = ir v (S(a)) st , 7r. v (o) = tt v (^(a)) 8 ', (A.7) 

respectively. Here st is the supertransposition operation defined by 

(A ab ) st = (-l^KM+M)^. (A.8) 

Note that in general ((A ab ) st ) st = (-l)^+^A ab ^ A afc . Let ist be the inverse operation of st such that ((A*&) s *) is * = 
((A a6 )"') st = A ab . Then 

(A ab y st = A ba = (-l)M+[ fc ] (,4 a6 ) st , (A.9) 

or A lst — 8A st 8, where 9 is a diagonal matrix with elements 8 ab = {— l)^S ab . 
By arguments similar to those in (20) , one can show that 

R v *' w {z) = {R vw {z)- 1 ) st \ R v -' w (z) = (R vw (z)- 1 Y t \ (A.10) 

It follows from the representations for R v ' W (z) and R V ' w (z) that for any pair of finite dimensional U q [Q^]- 
modules V and W, the R-matrix satisfies the following crossing-unitarity relations 

(((fl™^*)- 1 )'* 1 )- 1 )"* 1 = (Kv(q- 2p ) ® lw)((i? yW/ (z< Z -^)) stl ) s * 1 (7rv( g 2 ") ® 

((^(z)- 1 )^)- 1 )^ = (ly ® 7T VV ( (7 2 ''))((i? V ' W 'M ff ))^) S * 2 (ly ® TT^fg- 2 ")). (A.ll) 

Note also that 

(M? ±2p ) ® 7r w (q ±2 P))R vw (z) = R vw (z)(Tr v (q ±2p ) ® ^(9 ±2p ))- (A.12) 
Let us remark that if one uses the opposite coproduct and antipode of U q [G^], 

A(hi) =hi®l + l®hi, A(d) = d®l + l®<i, 
A(e t ) = e, ® + ® e u A(/i) = f l ® + q~^ ® /*, 

S(a) = -qPaq-P, Va = d, hi, <*, f t , (A.13) 

and denote the corresponding R-matrix by R(z), then the similar arguments as above give rise to the following 
crossing-untarity relations: 

(((R vw {z)- 1 )' tl )- 1 )" tl = (Ma 2 ') ® iw)(( J R VHr (««* ff ))* tl )" tl (Tv(?~ 2 ") ® V), 

(((^(z)" 1 ) 8 * 2 )- 1 )^ = (ly ® 7r w (g- 2 "))(( J R yw (zg-t«)) st2 ) st2 (lv ® tthV)). (A.14) 
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